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Abstract. Vertically aligned mono-domain nematic liquid crystal elastomers contract when heated. If a 
temperature gradient is applied across the width of such a cantilever, inhomogeneous strain distribution 
leads to bending motion. We modelled the kinetics of thermally-induced bending in the limit of a long 
thin strip and the predicted time- variation of curvature agreed quantitatively with experimental data from 
samples with a range of critical indices and nematic-isotropic transition temperatures. We also deduced a 
value for the thermal diffusion coefficient of the elastomer. 
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1 Introduction 

Liquid crystal elastomers (LCE) combine the long range 
orientational correlation of liquid crystals and entropically- 
driven polymer elasticity to give a range of exotic prop- 
erties such as spontaneous reversible shape changes and 
"soft elasticity" - deformation with no or little energy 
cost, leading to a variety of director instabilities under 
transverse extensions pQ. 

Nematic LCE possess the simplest uniaxial orienta- 
tional order and can be synthesized by incorporating rod- 
like anisotropic mesogenic groups into the strands of weakly 
cross-linked polymer networks. The order is character- 
ized by its principal axis, the nematic director n, and the 
scalar order parameter Q = (P^cosO)), which measures 
the mean orientation of mesogenic groups with respect to 
■ the director. Such an internal degree of freedom coupled to 
' the elastic body constitutes what is known as the Cosserat 
medium: the relative movement of crosslinking points pro- 
vides elastic strains and forces, while the director rotation 
causes local torques and couple-stresses - both intricately 
connected in the overall macroscopic response of the body. 
In fact, the physics of LCE is much richer than of notional 
Cosserat solids because (again due to the entropic nature 
of long polymer chains connecting the crosslinking points) 
rubbers are capable of very large shear deformations (be- 
ing at the same time essentially incompressible). Hence, 
one expects a variety of unique physical properties, espe- 
cially in the region of large deformations. However, in this 
work we shall explore only small local strains. 

Due to the coupling to the elastic body, the change in 
the degree of alignment of mesogenic rods leads to spon- 
taneous elongation or contraction of the whole network 
along n as constituent polymer chains become on aver- 
age more or less anisotropic (prolate or oblate depend- 



ing on the system). This direct coupling between physi- 
cal conformation and order parameter has been theoret- 
ically predicted a long time ago [2], and then compre- 
hensively demonstrated by simultaneous measurements of 
length and order using diffraction techniques [SUMS] . 

Landau theory predicts a 1st order transition into the 
isotropic phase as the nematic LCE is heated above its 
nematic-isotropic transition temperature Tjv/. This is based 
on the quadrupolar symmetry of the second-rank tensor 
order parameter of the nematic, which does not distin- 
guish between rods pointing "upwards" and "downwards" . 
However, in the elastomer network that was crosslinked 
in the aligned director state in order to obtain a per- 
manent monodomain nematic texture 6J, one does not 
find a discontinuous jump in the order parameter. Instead, 
frozen-in uniaxial stress leads to the supercritical continu- 
ous change of Q(T) across the transition, and with it - the 
continuous equilibrium uniaxial deformation of the mon- 
odomain nematic LCE. Depending on the degree of in- 
duced anisotropy of polymer chains forming the network, 
the magnitude of this deformation can be as high as 500% 

Spontaneous shape changes can also be induced opti- 
cally. Photoelastomers doped with rod-like groups, such 
as azobcnzcnc derivatives, which undergo trans- cis iso- 
merization on absorption of UV photons [HUH] , or carbon 
nanotubes, which respond to IR light [10 , are found to 
contract when irradiated at suitable wavelengths since lo- 
cal order is disrupted by the kinked dopant groups. Due to 
the high stroke and the equilibrium (reversible) nature of 
induced deformations, this now becomes an active area of 
engineering micro-optical mechanical systems (MOMS). 

Inhomogeneous deformations are of special interest since 
they see potential applications in photo- and thermal ac- 
tuators, detectors and sensors, microrheological valves and 
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pumps, as well as structures which can respond to their 
neighboring environment. Non- uniform deformations oc- 
cur when a spatial stress distribution is induced inside 
an elastomer. This could be achieved by irradiation, or 
by application of temperature gradient across the sample. 
Mathematical models [TT] have been proposed to predict 
equilibrium curvatures in unilaterally illuminated photoe- 
lastomer cantilevers with exponential attenuation. How- 
ever, in contrast to uniaxial contractions along the ne- 
matic director that have been well documented, no quan- 
titative measurements have been made so far on bend- 
ing curvatures. At the same time, while qualitative exper- 
iments on optically-induced deformations have reported 
time-scales varying from < 100ms [12 to ~ 1 — 10s [13l 
114] depending on incident intensity of light sources, the 
kinetic aspects of the bending motion have not been ad- 
dressed theoretically. 

This paper presents the first quantitative experimental 
study of the dynamics and kinetics of thermally induced 
bending in a nematic elastomer cantilever. We apply ra- 
diative heating to one side of cantilevers made from well- 
aligned monodomain polysiloxane side-chain elastomers, 
and measure the amplitude and time evolution of the in- 
duced curvature. We also develop a theoretical model, 
which predicts the reduced curvature of the cantilever as a 
function of time for cantilevers with with different critical 
exponents, transition temperature and maximum strain. 
A value for the thermal diffusion coefficient of the elas- 
tomer is estimated from matching the model predictions 
to the experiment. 



2 Experimental section 

Materials. All side-chain siloxane liquid crystalline elas- 
tomers, as well as their starting materials, were prepared 
in the Cavendish Laboratory following the procedures of 
Finkelmann et al. [6"lll5|. The polymer backbone was a 
poly-dimethylhydrosiloxane with approximately 60 Si-H 
units per chain, obtained from ACROS Chemicals. The 
pendant mesogenic group in sample A (NE-A) was purely 
4-methoxyphenyl-4-(l-buteneoxy) benzoate (MBB), while 
sample B (NE-B) contained of 70mol% of MBB and 20mol% 
of 4-alkeneoxy-4'-cyanobiphenyl (ACB), as illustrated in 
Fig. [TJ All networks were chemically crosslinked via the 
same hydrosilation reaction in the presence of commercial 
platinum catalyst COD, obtained from Wacker Chemie, 
with di- functional crosslinking group 1,4 di(ll-undecene) 
benzene (HUB) also synthesized in-house. In all cases the 
crosslinking density was 10 mol% of the reacting bonds 
in the siloxane backbone, so that on average each chain 
has 9 mesogenic groups between crosslinking sites. These 
materials are very well studied over the years; both have 
a glass transition around 0C and nematic-isotropic tran- 
sitions: T NI ks 87C for NE-A, and T NI w 101C for NE-B. 

Monodomain alignment. Mono-domain, aligned sam- 
ples of nematic elastomers were made by following the 
classical two-step crosslinking approach of Finkelmann et 
al. [6]. First we prepare partially crosslinked films in a 
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Fig. 1. Schematic illustration of the materials used in this 
work. Siloxane backbone chain with Si-H groups was reacting 
with 90 mol% mesogenic side groups and 10 mol% of flexi- 
ble crosslinking groups (HUB). Two materials differed in the 
composition of mesogenic groups: NE-A had 90mol% of MBB, 
while NE-B had 70mol% of MBB and 20mol% of ACB. 



centrifuge, highly swollen in toluene (2-3ml per lg of ma- 
terial), reacting for 25-35 minutes at ~75C before evapo- 
rating the solvent and suspending the samples under load 
in an oven for more than 5 hours at 120C to complete 
the second-stage crosslinking reaction. A careful study 
of reaction kinetics ensured that approximately 50% of 
crosslinks were established in the first stage of this prepa- 
ration. When a uniaxial stress is applied to such a partially 
crosslinked network, the uniaxially aligned state in the re- 
sulting nematic elastomer is established with the direc- 
tor along the stress axis. This orientation is then fixed by 
the subsequent second-stage reaction, when the remaining 
crosslinks are established. 

Following the original ideas of [B] and the present un- 
derstanding of the nature of polydomain nematic LCE 
16J, in all cases we performed the second stage crosslink- 
ing in the high-temperature isotropic phase: only in this 
way a good alignment and mechanical softness are achieved 
(in contrast to crosslinking in a stretched polydomain ne- 
matic phase, which results in topological defects and lo- 
calized domain walls frozen in the material). 

The mechanical history of the samples was eliminated 
by annealing in the isotropic phase for > 2 hours (~ 
130°C) followed by slow cooling. Precise measurement of 
variation in natural length L(T) with temperature was 
then made with a travelling microscope, which followed 
the end points of a sample that was suspended without 
load and heated at a slow rate of 0.33C/min in an insu- 
lated glass-front oven. 

Measurements of natural length L(T) variation with 
temperature, Fig. [21 were fitted to a model function L/Lq = 
1 + (3(1 — T/Tjv/) a , where L a is the constant length of 
samples in the isotropic phase. Obviously, such a superfi- 
cially critical behavior cannot be matched to experiment 
at the transition point itself (where supercritical effects 
take over), but it provides a very good continuous inter- 
polation of the data in the nematic phase. Fitting to the 
data gives a = 0.25, (3 = 0.843 and T NI = 359. 6K for NE- 
A, and a = 0.21, (3 = 0.3 and T m = 373.8K for NE-B. 

The thickness w of the two samples (which we shall re- 
quire in the cantilever analysis) was 0.351mm (NE-A) and 
0.368mm (NE-B) at room temperature. Elastomer sam- 
ples were cut into thin strips of approximate dimensions 
L x W x w — 5mm x 1mm xw (with differing thick- 
ness) and had one end vertically attached to a stand on an 
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Fig. 2. Curves of equilibrium uniaxial contraction of mon- 
odomain nematic LCE. The two samples (labelled on the plot) 
have different transition temperatures Tni and different chain 
anisotropy, leading to the 20% and 55% contraction respec- 
tively. Solid lines drawn through the data below Tni represent 
the analytical fit functions described in the text. 



adjustable platform. Sideways images of the strips (can- 
tilevers) were taken by a Sanyo VCB-3512T monochrome 
CCD camera (f = 9mm) with direct back lighting and digi- 
tally captured using software FTA32 by First Ten Angstroms 
Inc at a frame rate of 15fps. An Antex CS 16 W soldering 
iron provided heating. This soldering iron had a flat tip 
(cylinder of 4.5mm-diameter) which was providing uni- 
form radiative heating over the whole cantilever. Imaging 
of this tip also acted as a scale for confirming the thickness 
w of the samples by comparing dimensions on-screen. The 
soldering iron, which was allowed to equilibrate for 15 min- 
utes before each experiment, was mounted horizontally on 
a movable stand which can be slid to the desired posi- 
tion ( 2mm) in front of the mounted sample in under 0.2s, 
which marked the start of each kinetic measurement. Tem- 
peratures at the front and the back of the samples were 
measured with a thermocouple, however, not during the 
cantilever-bending experiment (but in a separate event of 
heating in exactly the same conditions). 

Movies of the bending motion were taken and each 
frame was analyzed both manually and using a MATLAB 
image-processing algorithm. Manually, the radius of cur- 
vature was obtained by superposing circles of various sizes 
on the outline of the curved sample using graphics soft- 
ware CorelDraw and adjusting until the circle of best fit 
was found, see Fig. [3] The automated MATLAB algorithm 
extracted the position of points along the curved edge of 
the sample in the image and fitted the set of points to 
the equation of a circle with variable radius R by a least- 
squares method; the optimized value of R was then out- 
put as the radius of curvature. In the end, our procedure 
was to analyze all images in an automated way, but then 
re-examine every anomalous point manually (because we 
found that our algorithm was not coping well with the 
cases of non- uniform curvature along the cantilever) . The 
outputs of this analysis were the values of normalized cur- 
vature w/R against time for each bending experiment. 



Fig. 3. (a) A composite image showing the extent of cantilever 
bending. On the left one can see the tip of soldering iron, (b) 
The scheme of manual analysis of cantilever curvature. 



3 Theoretical model 

Heat diffusion across a flat strip 

The problem of thermal diffusion in a flat sample ex- 
posed to a constant heat flux from one side is certainly a 
classical one. We give its brief account here in order to ex- 
pose the key parameters of the problem, required for the 
subsequent description of cantilever bending. Consider a 
ID diffusion of the scalar temperature field T(x, t) across 
the thickness of the cantilever, d{T — Dd 2 T, where the 
diffusion coefficient D = n/C is the ratio of thermal con- 
ductivity k (given by the heat flux definition J = —k8 x T) 
to specific heat capacity per unit volume C. We assume 
that the soldering iron acts as a source of constant flux J 
at x = while losses on the front (x = 0) and back (x — w) 
surfaces are taken to be proportional to the difference be- 
tween the temperature on the sides, Tf and Tb, and the 
ambient surrounding, To, i.e. —nd x T = J — j(Tf — To) at 
x = 0, and —nd x T — 7(T, — To) at x = w, taking flux to 
be positive to the right. 

Since the diffusion equation consists only of derivatives 
of the temperature field T and the boundary conditions 
are only sensitive to temperature differences, we can ho- 
mogenize the problem by considering the function 9 = 
(T — To) instead. Introducing natural variables x = x / w 
and r = Dt/w 2 , the problem can be recast as: 



d T e = d x e , 



with < d Y 



d x 9 = A(9 - e c ) at x = o 



-A9 atx = 
at r = 



(1) 



J/7 are the two essential 



where A = w^/k and C 
parameters of the problem. 

In the steady state d T 9 = 0. It being a 1-D problem, 
only a linear solution 9 — Ax + B could satisfy d 2 9 = 0. 
Letting the steady state front and back temperatures be 
and T f * and T^ respectively (equivalently 0j* and 6^), the 
steady-state temperature profile across the sample is given 
by 9 S = 9f - (9f - 9^)x- Parameters A and O c can then 
be expressed as 



A 



1) and c =(0f 



(2) 



which consist only of explicitly measurable quantities Tf , 
Td and T). Note that for a thin enough sample one expects 
The full time-dependent solution 



T* 



T b * and so A < 1. 



4 K. K. Hon, D. Corbett and E. M. Terentjev: Thermal diffusion and bending kinetics in nematic elastomer cantilever 




cantilever thickness. Since the length of a mono-domain 
nematic elastomer below T^i is locally given by L/Lq = 
1+/3(1 — T/Tf4i) a , the local strain distribution due to con- 
traction along the z-dircction (vertical in Fig. [3]) can be 
calculated. Taking the zero-strain state at ambient tem- 
perature Tq, we obtain: 



0.4 0.6 
Reduced depth x/w 

Fig. 4. Plots of 9 = T — To across the cantilever thickness x/w. 
Increasing curves are for the times Dt/w 2 — 0.1, 1, 3 and 100. 



can be obtained by superposing a series of different time- 
decay modes 9 = Ak sin(kx + 4>) exp(— k 2 r) on the steady 
state solution. Quantization conditions for k and <fi are 
determined by the initial/boundary conditions, giving the 
transcendental equations 

k n = At&n(nir/2 — k n /2) and tfi n = (nn/2 — k n /2) (3) 

for integer n. Fourier analysis of orthogonal modes in 0(x, t) 
gives the expression for coefficients 

An = i _ 1){n+m 2e c Mk /2) foroddn; (4) 

K n + sm(K n J 

_ / 1 ^«/2+i 26> c A k n cos(k n /2) - 2sin(fc»/2) 



kn 2- 



sin(fc n ) 



for even n. The full solution for the temperature across 
the sample of thickness w is therefore given by 



TCx,t) = T { * - (Tf-Tu)- 

w 



(5) 



^2 A n si 



sin k 



n=l 



There are two features of this solution that we need for 
our main problem. First of all, there is a characteristic 
time scale in the problem, given by the ratio w 2 /D which 
will allow us estimate the thermal diffusion constant in 
nematic LCE. Note that the (n = 1) mode in Eq.© has 
k\ ~ irA at A <C 1 and, therefore, this is the slow-decaying 
mode. Other modes have k n of order (n — l)w and decay 
fast, in practice, within a few seconds in our experiments. 

The second aspect of the solution T(x, t) is the rather 
smooth variation across the cantilever thickness. Figure 2] 
demonstrates the ^-dependence at different times, which 
justifies an essential simplifying approximation made in 
the next section, taking T(x) to be a simple linear func- 
tion connecting the two values Tf and 

Kinetics of cantilever bending 

Consider now a long thin strip of elastomer, preferen- 
tially contracted at the front and bent accordingly, due to 
unilateral heating from the side x = starting at t = 0. 
A temperature distribution T(x, t) is set up across the 



e{x) 



L{x) 
L(T ) 



(G) 







[1 + 13(1 - T IT NI Y] 



1 - 



T 
Tni 



- 1 



where L(Tq) is the sample length at ambient temperature. 
Assuming Young's modulus to be constant over the tem- 
perature and strain range of the experiment, this strain 
can be directly converted to stress a = Ee(x), at each 
depth x into the sample. As the elastomer deforms in- 
compressibly, a contraction in the z-direction would lead 
to transverse expansions along x and y. However, for the 
case of a long thin strip in which L 3> W > w, x- and y- 
curvatures can be safely ignored. 

Let x n be the position of a neutral plane [TT] . Mechan- 
ical equilibrium requires, in the absence of external forces 
and torques, that force and moment vanish across every 
cross-section of the cantilever. This means that all bending 
stress, E(x — x n )/R where R is the radius of curvature of 
the beam, is provided by the excess stress E[e(x) — e(x n )] 
with respect to this neutral plane. Two conditions repre- 
senting the balance of forces and torques are, as in pj] , 



WdxE 

R 



WdxxE — — 
R 



WdxE[e(x) -e{x n )\ (7) 
W dx x E[e(x) — e(x n )] 



where W is the width of the cantilever (the ^-dimension), 
Ww being the cross-section area of the beam. 

Assuming a linear temperature drop across the thick- 
ness, T(x) — Tf — (Tf — Tb)x/w, the spatial integrals can 
then be conveniently converted to over temperature via 
dT = — (Tf — T\,)dx/w to obtain 



(8) 



1 




/■Tf 


e{T)dT 

Tf - T b £[Xn} 


R 


(f ~ Xn ) 






"W x n \ 


-f 


(T f - T)e(T) dT e(x n ) 


-n( 


~ ~2J 




(Tf - T b ) 2 2 



Note that both the lateral width W and the Young mod- 
ulus E scale out of these mechanical balance equations. 
On elimination of x n and e(x n ) from these two equations, 
one obtains 



W 
R 



The problem could be solved completely if we had an an- 
alytical expression for s(T). Unfortunately, our interpo- 
lation formula is only applicable below Tjv/, while above 
the transition e is constant. The break at Tux calls for 
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a different mathematical treatment for the following two 4 Analysis of bending kinetics 



regimes: 

(i) When the temperature is below Tni everywhere in 
the sample, Tni > Tf > Tb, evaluating the integrals gives 
the dimensionless curvature: 



w 
R 



60 



1 



[1 + (3(1 - T /T NI ) a ] 2 + 3a + 



(10) 



T 



NI 



Tni 



+ 



1 - 



Tf-T b 



2±X 

Tni ) 
Tni) 



} Tni — Tf 
J T - T b 
T N i - T 
Tf-Tb 



Here Tf and T> are independently measurable time-depende 
functions. At all temperatures within this regime the cur- 
vature w/ R is a monotonically increasing function of time 
with the positive second derivative (convex function). 

(ii) When the Tf exceeds Tni, so that the phase transi- 
tion front is inside the sample, Tf > Tni > Tb, we obtain: 



w 
R 



6/3 



Tni — T> 



(3(1 - T /T N i) a } T-Tb 

T b \ n o + 2(r f -r J v J )/(T f -Tb) 
Tni) 2 + 3a + a 2 



(11) 



This expression represents a non-monotonic function of 
time, with the negative second derivative (concave func- 
tion) and the maximum curvature w/R followed by a de- 
crease when the majority of the sample becomes isotropic. 

The regime (ii) comes to an end when the temperature 
at the back of the cantilever, Tb reaches the transition 
point Tni, i- c - ah of the sample becomes homogeneously 
isotropic. At this point, evidently, w/R = 0. Calculat- 
ing the integrals in the regime when all of the sample is 
isotropic, Tf > T> > Tni, confirms that w/R = at all 
times. 

In our experiments the powerful heating flux has en- 
sured that the temperature rise was high so the interest- 
ing regimes (i) and (ii) occurred at relatively short times, 
when both Tf and Tb were well-approximated by a single 
exponential: 



T f 
T b 



T/ 



(T b * 



T )e- t / Tf 
T )e- t / Tb 



(12) 



where Tf and Tb are the effective thermal diffusion times 
at the front and back, respectively. It is expected that Tf 
and Tb would take similar but not identical values, since 
inspection of the full series for T(x, t) shows that the spa- 
tial coefficient of sin[k n (x/w — 1/2) + mr/2] is identical 
at x = and x = w for odd modes, but swaps sign for 
even n. In this single-exponential approximation, regime 
changes occur at t x = Tf ln(T f * - T )/(T f * - T NI ), when 
the front of the elastomer enters the isotropic phase and 
the curve w/R(t) has an inflection point; and at <2 = 
Tbln^b - Tq)/(T£ - T N i), when the elastomer becomes 
uniformly isotropic and hence returns to a state of zero 
curvature, with also a zero tangent. The values and gra- 
dients of reduced curvature match on both sides of t\ and 
ti, as is physically required. 



In both samples, the reduced curvature w/R variation is 
characterized by three distinct regimes. Immediately after 
the start of the experiment (heat flux on) there is a slow 
initial increase of curvature with time; this is followed by 
a sharp peak, after which the curvature rapidly dropped 
to zero. This trend corresponds to the regimes (i) and 
(ii) described in the theoretical model, but also could be 
qualitatively understood by considering the shape of the 
contraction curve of the nematic elastomer, L/Lq against 
T (see Fig. [5]). Initially, when the sample is relatively far 
away from Tni, the mechanical response to temperature 
n1 change is flat. The local strains induced at the front and 
back surfaces are small and similar in magnitude. How- 
ever, as the temperature approaches Tni, which happens 
first on the front of the cantilever, the internal stress gradi- 
ent is amplified by the increasing steepness of L(T) curve. 
This leads to cantilever curvature rising at an increasing 
rate. Curvature will be decreasing again after the front 
portion of the cantilever turns isotropic and stops con- 
tracting while the on-going contraction at the back re- 
duces the stress gradient. Finally, the sample returns to a 
stationary unbent state as it becomes uniformly isotropic. 

Figures [5] and [6] show the measured values of curvature 
w /R. The solid line in each plot is the fit by the theoretical 
model, Eqs. (fT0|) and (fTT|) . Evidently, the agreement is very 
good, both qualitatively and quantitatively - except in the 
final stages of sample un-bending, where the discrepancies 
are significant. This, however, should be expected because 
of the following two factors, one practical, the other to do 
with data analysis. 

In our model, we have ignored the fact that, as cur- 
vature increases, the lower end of the cantilever lifts and 
becomes closer to the heat source. Figure |3] shows this 
very clearly. This introduces a significant deviation from 
the theoretical assumption in the model, that the heat flux 
J is constant. It is clear that, after the point of maximal 
bending is reached, the real heat flux on the sample is 
inhomogeneous along the length of the cantilever (z- lo- 
cal coordinate). As a result, the far end of the cantilever 
(the part closest to the heat source) will become homo- 
geneously isotropic much earlier than the simple ID dif- 
fusion model would predict. In practice we see this very 
clearly, as in the final stages of heating cycle the far end 
of the cantilever is already straight, while the middle and 
near parts still have curvature remaining. Such buckling 
occurred more significantly in NE-B than in NE-A, as seen 
from the curvature plots, however in all cases both man- 
ual and algorithmic fitting was difficult and ambiguous as 
curvatures are no longer constant. 

The other discrepancy factor is inherent in the model 
itself. In order to obtain closed- form expressions for w /R(t) , 
Eqs. fTOf and (fTTj) , we had to use an interpolating formula 
for the underlying thermal contraction L(T). As it is clear 
from Fig. [21 this formula does not work in the immediate 
vicinity of the transition Tni, where all materials show 
differing degree of diffuse supercritical behavior. So it is 
not surprising that the model expressions deviate from 
the actual data when the temperatures of the sample are 
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Fig. 5. Experimentally measured values of reduced curvature 
w/R, for the NE-A cantilever, against time after the start of 
heating. The solid line is the fit by theoretical equations (|10l) 
and pip, with arrows showing where the regimes (i) and (ii) 
end. 
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Fig. 6. Experimentally measured values of reduced curvature 
w/R, for the NE-B cantilever, against time after the start of 
heating. The solid line is the fit by theoretical equations (|10l) 
and pip, with arrows showing where the regimes (i) and (ii) 
end. 



around, or slightly above, the notional Tjvj. These are the 
regions past the curvature peak, where the model "as- 
sumes" the local regions with T > Tni are fully isotropic, 
while in practice we know the contraction continues for 
2-3 more degrees. 

Among other ignored effects, which might become ap- 
parent at high temperatures (later times of the bending 
cycle) there is heating of surrounding air and lateral heat 
loss, which would require ambient temperature To to be 
time-dependent and a full 3-D treatment of heat diffusion 
respectively. All these factors could be eliminated or much 
reduced in impact. We could (and indeed have in some ex- 
periments) mount the heat source at an angle to minimize 
the effect of one cantilever end approaching it too close. 
We could also write a much more elaborate interpolation 
formula to account for the full continuous L{T) variation, 
and then proceed to calculate all integrals numerically. 
However, on reflection we have decided that the benefits 
of such improvements would not be worth the price of los- 
ing the simplicity of experiment and the ease of analysis. 



After all, the agreement of the model with experiment in 
the early regions of the bending process is excellent, as we 
expect when the material is in the nematic phase. 

The point of maximum curvature occurred at 4.3s for 
NE-A and at 3.3s for NE-B. Since we have measured the 
saturation temperatures T f * and T£ independently, as well 
as determined the parameters /3, a, Tjvj of the intrinsic 
thermal contraction curves for each material, the only 
two fitting parameters are the front- and back-relaxation 
times Tf and Tb, cf. Eq. (fl2|) . The best fits in Figs. [5] and 
|5] were achieved with Tf = 2.42s and Tb =2. 05s for NE- 
A, and with Tf =2. 13s and Tb —2. 5s for NE-B. It must 
be remarked that despite the experimental and theoreti- 
cal difficulties discussed above, the two-stage increase in 
the cantilever curvature shows good agreement with the 
experimental data. This proves that we understand the 
underlying physics of nematic elastomer cantilevers cor- 
rectly, and allows us to extract relevant material parame- 
ters from the fits. 

Parameters Tf and Tb can be viewed as a character- 
istic heating time for the sample. By comparison with 
the first terms of the thermal diffusion solution, Eq.([5]), 
they are expected to be of order w 2 /k 2 D. Steady state 
temperatures T f * and T£ were measured to be 375K and 
365K for NE-A, and 398K and 390K for NE-B, giving 
A = 0.89 and 0.133 respectively. The value of fei can 
then be obtained from solving the transcendental Eq.([3]) 
numerically (giving k\ — 0.419 for NE-A and 0.511 for 
NE-B). Rearranging, the thermal diffusion coefficient is 
given by D ~ w 2 /k\r. A typical value of D is therefore 
estimated to be ~ 1.5 • 10~'m 2 /s. We are not aware of 
any measurements of thermal diffusion in nematic LCE, 
but this estimate compares favorably with the literature 
values of D = 1.1 • 10~ 7 m 2 /s for a crosslinked silicone 
elastomer [17], a value also consistent with uncrosslinked 
silicone polymer melts [18] . 



5 Conclusions 

In this work we experimentally studied the time- variation 
of curvature of a long thin strip of aligned monodomain 
nematic elastomer, for two samples differing in their tran- 
sition (constitutive) behavior and shape dimensions. As- 
sociated theoretical analysis was able to quantitatively de- 
scribe the data and reflect all characteristic trends. Fitting 
the data allowed us to deduce characteristic time scales 
and estimate the thermal diffusion constant D of the silox- 
ane elastomer. 

To the best of our knowledge, this is the first quantita- 
tive study of thermal diffusion in nematic elastomers, as 
well as their cantilever bending due to induced inhomo- 
geneous strains arising from unilateral radiative heating. 
Such controlled and reproducible bending is an impor- 
tant physical effect underlying many engineering appli- 
cations. Perhaps more practically relevant is the photo- 
induced cantilever bending, where the local strains are in- 
duced due to photoisomerization reaction in azobenzene 
derivatives [TTl[T2lfTBlll4l[Trj] . Our experiments on photo- 
bending, analogous to the current work, will be reported 
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elsewhere. Nevertheless, thermal bending is a fundamen- 
tally important effect where one tests the details of con- 
tinuum mechanics, kinetics of local and global response, 
and the general understanding of nematic elastomer state. 
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